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Abstract 

We report on thermodynamic and optical measurements of the condensation process of ^He in 
two sihca aerogels of same porosity 95%, but different microstructures resulting from different 
pH during synthesis. For a base-catalyzed aerogel, the temperature dependence of the shape of 
adsorption isotherms and of the morphology of the condensation process show evidence of a disorder 
driven transition, in agreement with recent theoretical predictions. This transition is not observed 
for a neutral-catalyzed aerogel, which we interpret as due to a larger disorder in this case. 

PACS numbers: 64.60.-i, 64.70.Fx, 68.03. Cd, 67.70. +11 
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Condensation of fluids in mesoporous media, i.e. with pores sizes falling in the nm- 
/im range, has long been the subject of extended research, both fundamental and apphed 
From the point of view of statistical physics, this phenomenon offers an experimental 
realization of a phase transition in presence of disorder and confinement. From the practical 
point of view, it is widely used as a way to characterize the pore distribution. The most 
common case is that of sponge-like porous media, where the fluid is conflned in well-deflned 
cavities, either highly connected as, e.g., in Vycor glass, or with little or no connection, as 
in MCM41 silica. In the case where the dense phase of the fluid wets the substrate, one 
observes the progressive adsorption of a thin fllm at low vapor pressure, followed by a rather 
abrupt fllling at a pressure of order, but smaller, than the saturated vapor pressure. This 
phenomenon is coined capillary condensation, being driven by the surface energy between 
the liquid and the vapor. In the classical picture, the condensation pressure depends on 
this surface energy and the pore size through the so-called Kelvin equation [l| . The process 
is hysteretic, emptying occurring at a lower pressure than fllling. While this hysteresis is 
well understood at the level of the unique pore, its origin remains debated in the case of 
disordered porous media, with connected pores of different sizes jl, 2|. In particular, it has 
been recently shown by Kierlik et al [3J that it could result from the disordered nature of the 
porous media, rather than, as usually thought, of its geometry or topology ('pore-blocking'). 
In this new description, for a given pressure close to saturation, there exists a large number 
of metastable equilibrium states, the one being selected depending on the time history of 
the system, which gives rise to the observed hysteresis. 

As such, this description should apply to porous media of very different topology, silica 
aerogels ^. In these gels, the silica forms a complex arrangement of interconnected strands, 
resembling more a three-dimensional spider web than a sponge. Despite the positive cur- 
vature of the strands , which unfavors the fllm adsorption, aerogels flU below the saturated 
pressure, which, in a 'classical' picture, has been attributed to capillary condensation start- 
ing on the intersections between strands [sl. Experiments are however delicate to interpret 
as, due to their very low elastic modulus, the gels deform a lot upon adsorption of fluids 
with usual values of surface tension, such as nitrogen . This leads to the idea of studying 
these gels using Helium, which has a much lower surface tension. Such an experiment, albeit 
differently motivated, was flrst performed by Wong and Chan, about 15 years ago [?]. Their 
results showed capillary condensation, but, quite surprisingly, at a well deflned pressure 
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rather than on some range, as one would expect considering the wide distribution of scales 
in the silica aerogels. Together with specific heat experiments, this lead these authors to 
propose the existence of a genuine phase transition of helium inside aerogels, with a slightly 
depressed critical temperature with respect to the bulk case (Tc = 5.1957^). However, latter 
experiments by other groups performed on similar aerogels (95% porosity, gelation process 
with basic pH), did not confirm this conclusion. The sorption isotherms are hysteretic, and 



condensation occurs over a finite range of pressures [a, l9|, [id, lUj • Still, this range remains 
unexpectedly narrow, especially at low porosity [s], [l2| or low temperature Q. 



Both facts are in agreement with theoretical work [J, ll3| extending Kierlik et aFs de- 
scription to a 95% porosity aerogel, numerically synthesized using diffusion limited cluster 
aggregation (DLCA), a process known to account for the structure factor of base-synthesized 
aerogels. These studies show a disorder-induced transition similar to that predicted by the 
Random Field Ising Model at zero temperature At small porosity (or large disorder), 
filling takes place by a succession of small, microscopic avalanches, where a limited number 
of sites switch from gas to liquid at each pressure. The maximal size of these avalanches 
increases with the porosity, up to some critical porosity (corresponding to some critical 
disorder), where it diverges. Above this critical porosity, which increases with increasing 
temperature, filling involves a macroscopic avalanche, associated with a jump in the average 
density, at some well defined pressure. This implies a change of shape, from smooth to 
sharp, of the hysteresis loop, when the porosity is increased at a constant temperature or, 
alternatively, when the temperature is decreased below some critical value, at a constant 
porosity. In magnetic systems, an equivalent change of the hysteresis loop has been observed 
in Co/CoO disordered films 15|. Although the observed evolution of shape with temper- 
ature is an appealing argument in favor of the Kierlik et al vision of the condensation in 
disordered porous media, other explanations could exist as well. The aim of this paper is to 
further investigate the adequacy of the proposed mechanism, and its possible dependence 
on the aerogel structure, using combined optical and thermodynamic measurements. 



We used a one step process 



16] to synthesize two aerogels of same porosity (95%), but 
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different microstructures, by changing the pH from basic to neutral. Neutron 
jisl scattering experiments previously performed on similar samples reveal the sensitivity of 
the structure to the different kinetics of the sol-gel process in both cases. The size of the 
building silica units is larger for the base-catalyzed aerogel (BlOO, i.e. of density 100 kg/m^) 
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than for the neutral aerogel (N102), whereas the correlation length is larger for N102 than 
for BIOO. This results in a fractal structure over a much larger range for N102 than for 
BlOO, and, accordingly, a wider distribution of pores sizes. The correlation lengths of our 
own samples were deduced from their structure factor S{q) at zero q (as measured by light 
scattering) and their known porosities and fractal dimensions, using a theoretical expression 
of 5'(0) which has been experimentally checked 19(]. This gives values of order 10 nm for 
BIOO and 20 nm for N102. Sorption isotherms were performed on both samples at a number 
of temperatures, and, at the same time, the way they scattered light was studied, so as to 
gain information on the distribution of helium inside the aerogels, both at the microscopic 
and macroscopic levels. 

The experimental cell is a 4 mm thick copper block, traversed by a 20 mm diameter hole, 
and closed by two sapphire windows. Aerogels samples were sliced from 14 mm diameter 
cylinders. In the first experiment, we used a half disk of BIOO, about 2.7 mm thick. In 
the second experiment, we studied a 3.8 mm thick full disk of N102. The cell is mounted 
in a cryostat with 8 optical ports 45 ° apart, and is connected to the room temperature 
gas-handling system by two lines, one of which is used to fill or empty the cell through a 
regulated flowmeter, and the other is connected to a room temperature pressure sensor, with 
a resolution of 0.1 mbar. The flowrate is adjusted so that the stage of capillary condensation 
takes from 15 to 30 hours, which we checked to be long enough to get rate independent 
results. 

The helium mass in the aerogel is obtained by integrating the flowrate over time, and 
correcting for dead volumes as discussed in ref. 1 

phase ('liquid'), assuming that the vapor inside the aerogel has its bu 
density of the dense phase is that measured for the fully filled aerogel 

Figure [T](a) shows the isotherms obtained at four temperatures below the bulk critical 
temperature Tc for the two aerogel samples. The general trends are in qualitative agreement 
with the usual picture of capillary condensation; the hysteresis loop gets narrower and 
closer to the bulk saturation pressure Pgat as the surface tension decreases with increasing 
temperature. Also, it is smoother, and closer to Pgat for N102, consistent with the wider 
distribution of pores extending up to larger sizes. Quantitatively, we measure the position 
of the adsorption isotherm by APmax = P — Psat, at the point of largest slope of the 
adsorption branch. For a cylindrical pore of radius R, filling occurs at the stability limit 



It is converted into a fraction $ of dense 
k value and that the 
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of the adsorbed film, wliicli, for R much larger than the film thickness, is given by [20|, l21 | 
APmax = o'/ R-Pv/{pl — Pv) , where a is the helium surface tension, and pi and pv are 
the bulk liquid and vapor densities at saturation 22|]. For a collection of independent 
pores, one expects AP^ax to correspond to some typical value of the pores size. Although, 
a priori , aerogels cannot be so simply described, figure [2] shows that the equation above 
accounts for the temperature dependence of APmax for N102 in the whole temperature range, 
and for BlOO above 4.95 K, with characteristic 'pores' sizes i? ~ 35 nm and 20 nm. Both 
values are about twice the measured correlation lengths, a difference which is not surprizing, 
considering the complex microstructure of the samples. However, this simple approach fails 
to describe the position of the pressure plateau for BlOO at the two lower temperatures, which 
occurs at a smaller AP^ax than expected. Furthermore, in such an approach, the small width 
of the corresponding isotherms would correspond to a much narrower distribution of pores 
sizes than that could be deduced from the isotherms at 4.95 K and above. 

Optical measurements also show a markedly different behavior below 4.71 K, and above 
4.95 K, respectively. We illuminate the aerogel by a thin (100 pm wide), vertically polarized, 
He-Ne laser sheet under a 45° incidence with respect to its faces, and image it at 45°, 90°, 
and (for N102) 135° using CCD cameras. At any point of the intercept of the aerogel by 
the laser sheet, the brightness of the image is proportional to the intensity scattered by the 
aerogel in the direction of observation, and depends on the local strength and scale of the 
spatial fluctuations of the helium and silica densities. Figure [3] shows (a) images of BlOO 
for increasing helium fractions along the previous isotherms, and (b) the scattered intensity 
by three selected spots in the images. Quantitative analysis of these curves show that, 
below $ ~0.3, the scattered intensity is that expected for a thin film of helium covering 
uniformly the silica strands, whereas, beyond this value, it becomes much larger, implying 
the formation of liquid clusters correlated over distances larger than the correlation length 
of the silica HI]. The striking fact is that, while this increase of the scattered intensity occurs 
uniformly over the sample for the two larger temperatures, for the two lower ones, it is only so 
for fractions less than about 0.5. Beyond, liquid progressively invades the aerogel, resulting 
in the growth of dark regions. Comparison of figures [1] and [3]^b) shows that this invasion 
occurs along the vertical part of the isotherm. This could suggest that the bright and dark 
regions correspond to two coexisting phases at thermodynamic equilibrium. However, this 
is ruled out by the fact that the intensity scattered by the bright regions evolves with $, as 



5 



shown by figure [3]^b) at 4.71 K. On the other hand, the observations exhibit the behavior 
expected for macroscopic avalanches ; at any given spot, the transition to the dark state 
occurs discontinuously, when this spot is swept by the boundary between the bright and the 
dark regions. This is consistent with the occurrence of a macroscopic (but local) avalanche 
from a partly filled to a fully filled state. In this scenario, the macroscopic morphology, 
which is identical for 4.47 K and 4.71 K, and for different rates of condensation as well, must 
originate in macroscopic, heterogeneities of the average density of silica, which make some 
regions favor slightly more the liquid. This will make the pressure at which the avalanche 
takes place to vary continuously in space, so that the transition should occur successively 
at different locations in the aerogel, in qualitative agreement with the observations. 

In summary, for BlOO at 4.71 K and below, the isotherms exhibit a nearly vertical portion, 
their position is anomalous, and a significant fraction of the condensation takes place in an 
heterogeneous way. These facts are consistent with the occurrence of a disorder induced 
transition for BlOO between 4.95 K and 4.71 K. In contrast, for N102 at all temperatures, 
the isotherms are smooth, and the scattering is homogeneous, as shown in figure 111(a) : as 
a function of the condensed fraction $, the intensity scattered at 45° is the same along any 
vertical line of the picture, corresponding to a given depth, hence path length, inside the 
sample. There is a difference between points at different depths, but this (as well as the 
bump in the signal at $ ^ 0.5) can be explained by multiple scattering effects, which are 
stronger for N102 than for BlOO, due to the larger thickness and silica correlation length. 
Hence, no transition is observed for N102 in the temperature range where it is for BlOO. In 
the disorder induced transition scenario, this implies that the disorder does not only depend 
on porosity and is larger for N102 than for BlOO. We speculate that this is due to the wider 
distribution of pores sizes in this case. It would be interesting to see whether calculations 
performed on 95% aerogels numerically synthesized using the Diffusion Limited Aggregation 
algorithm appropriate for neutral aerogels would indeed confirm the absence of transition in 
this case. 

Finally, as shown by figure H] (b), the intensity scattered at 135° close to the entrance of 
the laser sheet, once corrected by the contrast factor ( (An)^, where An is the difference in 
refractive index between the liquid and the vapor) behaves roughly the same, whatever the 
temperature. This suggests that, at a given $, the distribution of liquid does not depend 
on temperature. This is not the case for BlOO below 4.95 K. As shown by figures [3](a) and 
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(b), the scattered intensity, at a given is smaller at 4.46 K than at 4.71 K, although the 
contrast is larger. This implies that the liquid must be distributed in smaller clusters. Here 
also, theoretical calculations would be needed in order to see whether such a temperature 
dependence follows from the disorder induced transition scenario. 

In conclusion, our results provide a strong evidence for the occurrence of a disorder in- 
duced transition in the condensation of ^He in base-catalyzed aerogels. Planned experiments 
using lighter aerogels should allow to check whether the transition moves to larger temper- 
atures, as the theory would predict. Here again, the very low index of refraction of helium, 
which allows quantitative single scattering measurements to be performed over a much wider 
range of conditions than for any other simple fluid, will be a unique advantage. This abil- 
ity to perform optical studies in transparent porous media could be used to study other 
problems as well, such as the origin of the hysteresis between adsorption and desorption. 

We acknowledge T. Herman, J. Beamish, F. Detcheverry, E. Kierhk, M.L. Rosinberg, and 
J. Phalippou for useful suggestions or discussions. 
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FIG. 1: (a) Isotherms at four temperatures, 4.47, 4.71, 4.96 and 5.08 K, for two aerogels of 
porosity 95%, BlOO (black) and N102 (grey). The hquid fraction <I> in aerogel is plotted versus 
AP = P — Psat] foi" BlOO, there is a clear change of shape with temperature. 
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FIG. 2: Position AP^ax of the point of maximal slope in the adsorption branch versus the temper- 
ature dependent ratio apv/{pL — Pv)- The straight lines correspond to pore diameters (see text) 
of 18 and 36 nm. 
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FIG. 3: Images of BlOO observed at 45° from the incident laser sheet for helium fractions of 
0.1,0.5,0.7,0.9, and 1.0. The curves give the scattered intensity, referred to the empty aerogel 
situation, for the three regions labelled (a), (b), (c). 
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FIG. 4: Intensity scattered by N102 as a function of the condensed fraction ^ : (a) at 4.47 K 
and 45°. The signal, normahzed by the sihca contribution, corresponds to five regions at three 
different depths inside the aerogel disk. Differences between curves for different depths (and path 
lengths) are consistent with multiple scattering effects combined with an homogeneous distribution 
of the liquid clusters within the sample, in contrast to BlOO at the same temperature; (b) at 135°, 
close to the entrance of the laser sheet. The silica contribution Isi has been subtracted, and the 
result normalized by Isi and the optical contrast squared. The signal is similar at all temperatures, 
suggesting that the spatial organization of the liquid only depends on not on temperature. 
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